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Abstract 



We study the semileptonic decays B — > Xse^e~, B -^ XgfJ, ^i~ in generic super- 
symmetric extensions of the Standard Model. SUSY effects are parameterized using the 
mass insertion approximation formahsm and differences with the Constrained MSSM 
C^ ' results are pointed out. Constraints on SUSY contributions coming from other pro- 

cesses {e.g. b -^ 57) are taken into account. Chargino and gluino contributions to 
photon and Z-mediated decays are computed and non-perturbative corrections are 
considered. We find that the integrated branching ratios and the asymmetries can 
be strongly modified. Moreover, the behavior of the differential Forward-Backward 
asymmetry remarkably changes with respect to the Standard Model expectation. 



1 Introduction 

One of the features of a general low energy supersynimetric (SUSY) extension of the Standard 
Model (SM) is the presence of a huge number of new parameters. FCNC and CP violating 
phenomena constrain strongly a big part of the new parameter space. However there is still 
room for significant departures from the SM expectations in this interesting class of physical 
processes. It is worthwhile to check all these possibilities on the available data and on those 
processes that are going to be studied in the next future. In this way it is possible to indicate 
where new physics effects can be revealed as well as to establish criteria for model building. 

In this work we want to investigate the relevance of new physics effects in the semileptonic 
inclusive decay B -^ X^f+f". This decay is quite suppressed in the Standard Model; however, 
new 5-factories should reach the precision requested by the SM prediction [|l| and an estimate 
of all possible new contributions to this process is compelling. 

Semileptonic charmless B decays have been deeply studied. The dominant perturbative 
SM contribution has been evaluated in ref. ^ and later two loop QCD corrections have been 
provided [Q, Q. The contribution due to cc resonances to these results are included in the 
papers listed in ref. ||5|. Long distance corrections can have a different origin according to 
the value of the dilepton invariant mass one considers. 0{l/ml) corrections have been first 
calculated in ref. @] and recently corrected in ref. [^ ||. Near the peaks, non-perturbative 
contributions generated by cc resonances by means of resonance-exchange models have been 
provided in ref. [0, PJ. Far from the resonance region, instead, ref. |jTO| (see also ref. ||11|| ) 
estimate cc long-distance effects using a heavy quark expansion in inverse powers of the 
charm-quark mass {0{l/ml) corrections). 

An analysis of the SUSY contributions has been presented in refs. [0-[|r^ where the 
authors estimate the contribution of the Minimal Supersymmetric Standard Model (MSSM). 
They consider first a universal soft supersymmetry breaking sector at the Grand Unification 
scale (Constrained MSSM) and then partly relax this universality condition. In the latter 
case they find that there can be a substantial difference between the SM and the SUSY 
results in the Branching Ratios and in the forward-backward asymmetries. One of the 
reasons of this enhancement is that the Wilson coefficient C7{Mw) (see section ^ for a precise 
definition) can change sign with respect to the SM in some region of the parameter space 



while respecting constraints coming from b ^ S'j. The recent measurements of 6 ^ 57 [16 
have narrowed the window of the possible values of Ci{Mw) and in particular a sign change 
of this coefficient is no more allowed in the Constrained MSSM framework. Hence, on one 
hand it is worthwhile considering B -^ Xsi~^i~ in a more general SUSY framework then just 
the Constrained MSSM, and, on the other hand, the above mentioned new results prompt 



us to a reconsideration of the process. In reference |jT^ the possibility of new-physics effects 



coming from gluino-mediated FCNC is studied. Effects of SUSY phases in models with 



heavy first and second generation sfermions have been recently discussed in ref. |]18 . 

We consider all possible contributions to charmless semileptonic B decays coming from 
char gino- quark- squark and gluino-quark-squark interactions and we analyze both Z-boson 
and photon mediated decays. Contributions coming from penguin and box diagrams are 



taken into account; moreover, corrections to the mass insertion approximation (see below) 
results due to a light iji are considered. A direct comparison between the SUSY and the SM 
contributions to the Wilson coefficients is performed. Once the constraints on mass insertions 
are established, we find that in generic SUSY models there is still enough room in order to 
see large deviations from the SM expectations for branching ratios and asymmetries. For our 
final computation of physical observables we consider NLO QCD evolution of the coefficients 
and non-perturbative corrections (0(l/m^), 0(l/m^),...), each in its proper range of the 
dilepton invariant mass. 

Because of the presence of so many unknown parameters (in particular in the scalar 
mass matrices) which enter in a quite complicated way in the determination of the mass 
eigenstates and of the various mixing matrices it is very useful to adopt the so-called "Mass 
Insertion Approximation" (MIA) |]T9[ . In this framework one chooses a basis for fermion and 
sfermion states in which all the couplings of these particles to neutral gauginos are fiavor 
diagonal. Flavor changes in the squark sector are provided by the non-diagonality of the 
sfermion propagators. The pattern of fiavor change is then given by the ratios 

/r/x _ i^ij)AB /.N 

sq 

where {m(j)\^ are the off-diagonal elements of the f = u,d mass squared matrix that mixes 
fiavor i, j for both left- and right-handed scalars {A, B =Left, Right) and M^g is the average 



squark mass (see e.g. [^). The sfermion propagators are expanded in terms of the (5s and 
the contribution of the first two terms of this expansion are considered. The genuine SUSY 
contributions to the Wilson coefficients will be simply proportional to the various (5s and 
a keen analysis of the different Feynman diagrams involved will allow us to isolate the few 
insertions really relevant for a given process. In this way we see that only a small number of 
the new parameters is involved and a general SUSY analysis is made possible. The hypothesis 
regarding the smallness of the 5s and so the reliability of the approximation can then be 
checked a posteriori. 



Many of these (5s are strongly constrained by FCNC effects PO, pi], B^j or by vacuum 



stability arguments p3| . Nevertheless it may happen that such limits are not strong enough 



to prevent large contributions to some rare processes. For instance it has been recently found 
in ref. [^ that the off-diagonal squark mass matrix elements can enhance rare kaon decays 
by roughly an order of magnitude with respect to the SM result. 

The paper is organized as follows. In sect. ^ we define the operator basis, the basic for- 
mulae for the BR, the Forward-Backward asymmetry and the non-perturbative corrections. 
Sect. ^ and sect. ^ treat chargino and gluino contributions in the mass insertion approxima- 
tion. The light Ir corrections are presented in sect. ^ Constraints on (5s are discussed in 
sect. 1^ and final results and conclusions are drawn in sect. |^ and ^. 
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Figure 1: Some of the relevant penguin diagrams for b -^ si~^i~ . The bubble indicates a 
Mass Insertion. Diagrams A,B are based on chargino interaction. Diagrams C,D consider 
gluino interactions. 
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Figure 2: Relevant box diagram for b -^ si^l . The bubble indicates a Mass Insertion. 



2 Operator basis and general framework 

The effective Hamiltonian for the decay B -^ Xgi'^i' in the SM and in the MSSM is given 
by (neglecting the small contribution proportional to K*gKub) 
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q. This Hamiltonian is known at next-to- leading 



K is the CKM-matrix and qL{R) 

order both in the SM [^ ^ and in the MSSM |14, 15]. We find that the most general 
low-energy SUSY Hamiltonian also contains the operators 
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However it is shown in following sections that the contribution of these operators is neg- 
ligible and so they are not considered in the final discussion of physical quantities. SUSY 
contributions to other operators are negligible because they influence our observables at an 
higher perturbative order. 

With these definitions the differential branching ratio and the forward-backward asym- 
metry can be written as 
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where s = (pf+ +p£-)^/m^, 6 is the angle between the positively charged lepton and the B 
flight direction in the rest frame of the dilepton system, f{z) and k{z) are the phase space 
factor and the QCD correction factor {z = mc/mb) that enter r{B -^ XcCv) and can be 
found in refs. 0, ^. C^^^s) includes all the contributions of the operators Qi — Qq and Qg 
and its complete definition for the SM and MSSM can be found again in refs. P, §, |14| . 

In the literature the energy asymmetry is also considered |]n| but it is easy to show 
that these two kind of asymmetries are completely equivalent; in fact a configuration in 
the dilepton c.m.s. in which ^^ is scattered in the forward direction kinematically implies 
Ei+ < El- in the B rest frame (see for instance ref. 0). 

It is worth underlying that integrating the differential asymmetry given in eq. (^ we do 
not obtain the global Foward-Backward asymmetry which is by definition: 
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where £1, and £^ stand respectively for leptons scattered in the forward and backward 
direction. 



To this extent it is useful to introduce the following quantity 
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whose integrated value is given by eq. (^. 
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Eqns. d^) and (|^) have been corrected in order to include several non-perturbative effects. 
First of all 0(l/m^) effects have been estimated by 0, 
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where 



^U 



3 - 8^2 + 24/ - 24^^ + hz^ + 24^^ log [z] 



Ai and A2 are the two parameters that appear in the Heavy Quark Expansion (HQE). While 
the value of A2 is quite well-established, (A2 = [M'b* — ^b)/4), Ai is not yet well known. 
In ref. ^ Ai is estimated as Ai = -0.52 ± 0.12 GeV^ and in ref. [g7| Ai = -0.10 ± 0.05 
GeV^. In what follows we consider the weighted average of the two results Ai = —0.16 
GeV^. As was pointed out by the authors of ref. [^ these corrections are no more valid 
near the endpoint region, s — *> 1, where they diverge because of the breaking down of the 
HQE. Following some recent analyses we have stopped the BRs corrections given in eq. (^ 
at s = 0.78 (see ref. 0) and the ones in eq. (p!0|) at s = 0.4 (see ref. 0). 

In order to account for the corrections to the parton model approximation in the high s re- 
gion ref. 10] and ref. ||^ adopt two different approaches. The former considers a Fermi-motion 
model and the latter invokes the Heavy Hadron Chiral Perturbation Theory (HHChPT). A 
discussion about the usefulness of the Fermi-motion model for semileptonic charmless B- 
decays is beyond the scope of this paper. In order to have a model independent description 
of the high energy region of the spectrum we have considered the HHChPT corrections. 

For 0.88 < s < s^ax = {Mb — MkY /ml = 0.99 the branching ratio is dominated by the 
exclusive decays B -^ Ki'^£~ and B — > K7r£~^i~; in ref. [Q is shown that the contribution of 
the latter is completely negligible. In the following we report the expression of the branching 
ratio, valid in the interval of the spectrum given above, computed in the HHChPT framework 
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and the definitions of /i, f+,Cf^, v ■ px can be found in ref. [§]. Moreover, we have put 
g=0.5 according to the theoretical estimate given in ref. |2^. In the intermediate region 
0.78 < s < 0.88 we have interpolated the obtained results. The form factors /+ and a^ 
can be computed also using other methods ( QCD Sum Rules ||29|, Ligh Cone QCD Sum 
Rules |3^, QCD relativistic potential model [^) but the HHChPT approach is preferable 
in the endpoint region of the spectrum. 

The asymmetries receive no contribution from the single kaon mode B — » Ki~^i~ and the 
endpoint of their spectrum is fixed, instead, by S ^ K7i£~^i~ at s = {MB — MK — MT,y/ml = 
0.93. In the region 0.7 < s < 0.93 we use the parton model result because these asymmetries 
have not been computed yet in the HHChPT framework. 

Finally, the following 0{l/ml) corrections occur for s < 0.2 (see e.g. HIT 
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The 0{l/ml) and 0(l/m^) corrections to Aps can be easily computed because of the 



following relation 
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All the effects coming from the mass insertion approximation can be included in formu- 
lae (|§-|T^) writing the coefficients Cj, Cg^{s), Ciq as 



Cr 



^SAI 



CI 



eff/ 



(Cfis))^'' + (Cf )^'"^ + {CfY", 



c 



10 



'-'10 + 



^Diag . ^MI 



ao 



'10 



(16) 



Diag 



where all the contributions are evaluated at the Mb scale and the various Q summarize 
all the contributions coming from graphs including SUSY Higgs bosons and sparticles in 
the limit in which we neglect all the mass insertion contributions (they would be the only 
SUSY diagrams if the scalar mass matrices were diagonalized by the same rotations as those 
needed by the fermions). The explicit expressions for Q ^"'^ can be found in ref. [P |. 

The Feynman diagrams with MI relevant for b -^ sl'^Cr are drawn in figs |l|-0. We have 
considered gluino-like and chargino-like contributions with both single and double mass 
insertions. 

Both photons and Z-bosons can mediate the decay. Usually one finds that Z-boson 
contributions are dominant in those graphs where an "explicit SU{2)l breaking" is provided. 
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Table 1: Central values of physical constants used in the phenomeno logical analysis 



i.e., both Left and Right squarks are present in the same loop. In the latter cases the photon 
cannot feel any gauge-symmetry-breaking and its contribution to the Wilson coefficients is 
suppressed by a factor m%/M]g with respect to the Z-boson one. For Mgg ~ 300GeV, this 
factor amounts to about an order of magnitude. On the other hand if the graph does not give 
any SU{2)l explicit breaking we are in the opposite situation and the Z-boson contribution 
is suppressed by a factor m1/m\ ~ 3 ■ 10~^. Moreover, a general feature of 7-mediated four 
fermion contributions is that, for high average squark masses, they decouple much faster 
than in the Z-boson case. This can be understood simply using dimensional arguments. 
While Wilson coefficients for Z-boson mediated four-fermion interactions are proportional 
to A/(m|Mg^ ), the same coefficients must be proportional to A/M^ for the 7 (A here is 



a generic off-diagonal element of the sparticle squared mass matrices and it cannot rise as 
fast as M^ for high values of M^q). Thus photon graphs can compete with Z-boson graphs 
if the sparticle spectrum is not too heavy. 

Finally the value of the physical constants we use is reported in table |I|. 

In the following subsections we describe in detail the contributions of each graph of 
figs. ^|. 



3 Chargino interactions 



In the weak eigenstates basis the chargino mass matrix is given by 



Mv 



M2 y2MH/sin/3 

\/2M\Y cos (3 jj, 



(17) 



where the index 1 of rows and columns refers to the the wino state, the index 2 to the 
higgsino one, /x is the Higgs quadratic coupling and M2 the soft SUSY breaking wino mass. 
In order to define the mass eigenstates the unitary matrices U and V which diagonalize M^ 
are introduced, 

diag(M^„M^J = t/*M^V^+. (18) 



After the rotation to mass eigenstates it is always possible to speak of wino-quark-squark or 
higgsino-quark-squark interactions. In order to identify the wino and higgsino states from 
the chargino ones it is sufficient to pick up the right elements from the U and V matrices. 
To be clear we write them explicitly for the cases of interest in the super-CKM basis. The 
wino-quark-squark, diWu, vertex is {d and u are a generic down-quark and up-squark) 

ig2 J2 ^kj J2 ValdjLXaUkL + h.c. , (19) 

j,k=l,2,3 a=l,2 

and the higgsino-quark-squark, dHu, vertex is 

^ E K*kj E [(^")fcK2 djLXaUkR - {\'^)kU*2 djRXaUkL] + h.C. (20) 

i,fc=l,2,3 a=l,2 

where K is the CKM-matrix, 02 = e/smOw, A'' = ^^" = ^^"^^ and A'^ = ^/2M1 

' ^sm/:/ y/2smpMw t^cos^ 

are the Yukawa matrices for the up- and down-quarks. 

Chargino graphs can contribute to the decay via both single and double insertions (see 
figs. |I]-0). The double insertion is particularly convenient if the corresponding 6s are not 
very constrained p3. In the following subsections we examine both the cases. In the case 



of a single insertion approximation both 7- and Z- mediated decays are considered. 
In all what follows our results for the integrals are written in terms of the functions 

To get a feeling with numbers it is sufficient to say that for a = b = 1, 

P,,fc(l,l) = /3s[l + z,l+j] (22) 

where Pe is the Euler /3-function. 



3.1 Single mass insertion — Z 

The Z-boson mediated decay can proceed in two ways depending on the type of chargino- 
quark-squark vertices we consider. If an explicit SU{2)l breaking on the squark line of 
fig. [l|A is required we must take both an higgsino and a wino vertex. In this way we get a 
contribution to the Wilson coefficients 



1 - 4 sin2 Ow ^ ^" ^"^--^^^2 Kt% 4 sin^ Ow ,,,tt,2 
I U*^Ujiy/XiXjPii2{xi, Xj) + Vi*^VjiPni{xi, Xj) - -6ijPo2iixi, Xj) [ , (23) 

where Xi = M^JM^ . This diagram, however, is exactly null in the limit in which U, V 
approximate the identity matrix and so it is negligible for high M2. 

9 



With two wino-quark-squark vertices we obtain 

i — 4sm t/vK -f^ts 4 sm (7^4/ j -^j^ 2 

{t/*if/ji^/x^Pii2(xi, Xj) + V*^jiPix\{,Xh ^j) - ^ijPm\{^i, ^j)\ (24) 

where we have retained only the contribution which arises because of the exphcit SU{2)l 
breaking (with a double wino-higgsino mixing in the wino line); in fact eq. (^) is null in 
the limit of diagonal chargino mass matrix. 

Graphs with two higgsino-quark-squark vertices are suppressed with respect to these ones 
by Yukawa or CKM factors. 

3.2 Single mass insertion — 7 

The contributions of the 7 penguin with two wino vertices are 

/If 2 1 L^* fi 

C7 = -(^23)lLJ^^J^^ E V^iV,\ |-P222(x.,X.) + Prs2ix.,X.) 
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sq 

The contributions of the 7 penguin with an higgsino and a wino vertices are 
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^^■'^sq 'J 92 J^ts ^b j=i_2 *- ^ ^ 

3.3 Single mass insertion — box 

Finally we compute the contributions which come from chargino box diagrams of fig. |^. 
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In the wino exchange case the result is 

K* M^ 1 
C, = -Cio = iS23)LL-^T^-^^ E iV,lV,^V,,V;,)fix,,x„x,) . (27) 

where 

1 /•! /•! /•! l/z(l-2)2 



^^ ^ ^ 2 7o io ^ Jo Ml 



z) + Xu{l - y){l - z) + z{xiX + Xj(l - x))] 

(28) 



and X£> = M?/M^g. 

If the wino-bottom-stop vertex is replaced by an Higgsino-bottom-stop one we obtain 

K* M^ A 

C9 = -C,o = -iS23)LR-^j^ . L E iV:,V,,VaV;,)f{x,,x,,x,). (29) 

3.4 Double mass insertion — Z 



It was recently pointed out |24| that a double mass insertion can provide a great enhancement 
of the SUSY contribution to the decay width, at least in the K-system case, if the Ss are not 
very constrained. 

For B decay we obtain contributions from this graph to Cg and Cio, 
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4 Gluino interactions 

The main contribution of this kind of interactions comes from the graphs drawn in fig. 01!, D. 
In what follows we analyze the single and double mass insertion cases. 

4.1 Single mass insertion — 7 

The corrections to the coefficients in the photon mediated decay case are: 
V2 1 N^ -1 Tia, 
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The term proportional to the gluino mass in eq. (pi]) seems to be strongly enhanced with 
respect to the others. However the mass insertion which enters the diagram is also strongly 
constrained from b ^ S'~f [^]. 

4.2 Single mass insertion — Z 

The only relevant contributions to the Z-boson mediated decay width come from diagrams 
in which the Z feels directly the breaking of SU{2)l. According to the argument of section § 
all the diagrams that do not respect this condition are suppressed with respect to the photon 
mediated ones and can be neglected. However, for penguins containing a gluino, an explicit 
SU{2)l breaking can be provided only with a double MI. If only one MI is considered, 
Z-mediated decays are completely negligible with respect to the 7-mediated ones. 

4.3 Double mass insertion — Z 

For completeness we report here also the result obtained performing a double mass insertion 

^- - ^ — 2Ari^ "^^"^"'"^ • ^''^ 

5 Light tR effects 

In the Mass Insertion Approximation framework we assume that all the diagonal entries of 
the scalar mass matrices are degenerate and that the off diagonal ones are sufficiently small. 
In this context we expect all the squark masses to lie in a small region around an average 
mass which we have chosen not smaller than 250 GeV. Actually there is the possibility for 
the tR to be much lighter; in fact the lower bound on its mass is about 70 GeV. For this 
reason it is natural to wonder how good is the MIA when a tn explicitly runs in a loop. 

The diagrams, among those we have computed, interested in this effect are the chargino 
penguins and box with the (^23)^^ insertion. To compute the light-t/? contribution we adopt 
the approach presented in ref . ||3^ . There the authors consider an expansion valid for unequal 



in the 


gluino 


penguin 


c. 






1- 


4sin^ 


6w 




1- 


4sin^ 


9w 



diagonal entries which gives exactly the MIA in the limit of complete degeneration. 

The new expressions for the contributions to the coefficients Cg and Ciq are the following. 
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Cliargino Z-penguin with both an higgsino and a wino vertex: 

- CiQ= i^23)LR-^-r^-r—r^ J2 ^ii^/2 ^ 



1 - 4 sin^ Ow " ""' 92 Kl 4 sin^ Ow ,,,tt,2 

where xr = {mi iMsqY and the functions j{x,y,z) and k{x,y,z) can be found in 



ref 32 



Chargino 7-penguin with both an higgsino and a wino vertex: 



M^2XtK. 



X 



tRJ 5 



where 



C9 = is-23)Ln^~-^Y.v^^y:2 97ix.., 

^sq ^92 J^ts i=l,2 

l\/r2 1 \ jy-* 

C, = {52,)LRj^~-^Y.y^^yr2 9M.Xi,) (34) 

^^^sq ^92 J^ts j=l,2 



9i{x,y) = — — 



frix) 



52 - 1532; + 144a;2 - 43x3 + (36 - 5Ax + 12a;3) log(a;) 



6(-l + x; 



4 



,,, -8 + 3x+l2x^ -7x^ + Qx(-3 + 2x)\og(x) 
^'^'^ = 6(-l+x)^ ■ ^''^ 

Chargino box with an higgsino vertex: 

C9 = -C^o = \iS-,)LR^^ .\, Y.iyayny^ly;2)k{x^,X„X„XiJ (36) 

4 J^ts ^^^sq fl'2Sin Uw jj=l,2 

where k{x,y,z,t) is defined in ref. ||3^ . 



All the above formulas reduce exactly to those presented in sect. ^ in the limit Xt^ = 1. 



6 Constraints on mass insertions 

In order to establish how large the SUSY contribution to 5 ^ Xgi^i^ can be, one can 
compare, coefficient per coefficient, the MI results with the SM ones taking into account 
possible constraints on the 6s coming from other processes. 

The most relevant 6s interested in the determination of the Wilson coefficients C7, Cg 
and Cio are [6^3)11, {62z)lr, ('^33)i?x, {^tzlLL and ((5|3)lr. 
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Vacuum stability arguments regarding the absence in the potential of color and charge 
breaking minima and of directions unbounded from below [^ give 



{51)lr < m 



2M| + 2Mf 

sq 



'M. 



(37) 



sq 



For Mgq < 300GeV this is not an effective constraint on the mass insertions. 

A constraint on {62^) ll can come from the possible measure of AMb^- 

In fact the gluino-box contribution to AMb^ p2[ is proportional to {62^) 
instance ref. |2^). A possible experimental determination of AMb^, say 



2 

23 J LL 



AMb,, < 30 ps"^ 



fsee for 



(38) 



would imply that 



Re((5|3)? , < 0.5 



'23JLL ^ "-^ (39) 

for squark masses about 250 GeV. Moreover the LL up- and down-squark soft breaking 
mass matrices are related by a Cabibbo-Kobayashi-Maskawa rotation 



[^sqlLL — -"^ V^-'sqlLL 



K\Ml?rrK 



SO that the limit (|39| ) would be valid for the up sector too: 



M823)Il < 0.5 . 



(40) 



(41) 



Some constraints come from the measure oi B ^ X^'-y. The branching ratio of this 
process depends almost completely on the Wilson coefficients C-j and C'-j which are 
proportional respectively to ((^23)^^2 or rl and ((523)ii- The most recent CLEO estimate 
of the branching ratio for i? ^ X^ 7 is [^ 

(42) 



Bexp(5 ^57) = (3.15 ± 0.35 ± 0.32 ± 0.26) ■ 10" 



where the first error is statistical, the second is systematic and the third comes from 
the model dependence of the signal. The limits given at 95% CL. are [Q: 



2.0 10-4 < Be..p(fi 



s 7) < 4.5 10" 



(43) 



We can define a Cf'{MB) as 



Cf{MB) 



BexpjB ^ g 7) 

Kcb J Trg{z) 



(44) 



where F can be found for instance in ref. |^. Considering the experimental limits 
given in eq. (|43|) we find 

0.28< |Cf (Mb)| <0.41. (45) 
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Actually |Cf (Mb)|2 = \Cr{MB)\'^ + |C7^(Mb)|2 and the constraint given in eq. (H) 
should be shared between the two coefficients. However in order to get the maximum 
SUSY contribution, we observe that in physical observables Cj does not interfere with 
Cy, the C7C9 term is suppressed by a factor mg/mb with respect to the C7C9 one and 
CjCq is numerically negligible (in fact Cg is much smaller than Cg). For these reasons 
we choose to fill the constraint of eq. (^5]) with Ct{Mb) alone. 

The bounds (^Sp are referred to the coefficient evaluated at the Mb scale while we are 
interested to the limits at the much higher matching scale. After the RG evolution 
has been performed we find that for an average squark mass lower than 1 TeV, the 
MIA contribution alone with a suitable choice of (5s, can always ffi the experimental 
constraints. 

Thus, since we are interested in computing the maximum enhancement (suppression) 
SUSY can provide, we can choose the total C^ {Mb) anywhere inside the allowed 
region given in eq. (P5| ) still remaining consistent with the MIA. 

The limit we get for ((523)l_r is of order 10~^ and this rules out Z-mediated gluino 
penguins contributions to Cg and Cio- 

For what concerns ((523 )ll ^^ ^"^^ ^^^^ ^^^ constraint changes significantly according 
to the sign of C^ {Mb)- In this case it is important to consider both the positive 
and negative region as this delta can give a non negligible contribution to Cg and 
CiQ. The limits depend on the choice of the parameters in the chargino sector; the 
numerical results given below are computed for Mgq ~ 250 GeV, /z ~ —160 GeV, 
Mj> ~ 50 GeV, tan/5 ~ 2 (in sect. ^ we will show that these are the conditions under 
which we find the best SUSY contributions). Considering the positive interval we find 
—0.7 < {623)11 < —0.5 while in the negative one \{S23)ll\ < 0.1. 

Finally a comment on the 6s coming in graphs with a double MI is in order. 

Given the constraints on {S23)lr one can see that the gluino-penguins with a double 
MI give negligible contributions to the final results even if {S^3)rl is of order 0{1). 

A {6^3) BL of order 0{1), can give rise to light or negative squark mass eigenstates. In 
particular a light ti would contribute too much to the pvK-parameter. Eventual model 
dependent cancelation can provide an escape to these constraints. In any case the nu- 
merical value of these contributions is not particularly important for the determination 
of physical observables. Since we want to provide a model independent analysis we 
prefer not to consider in our final computation these double insertion graphs and we 
present them only for completeness. 

Contributions with three mass insertions are suppressed due to small loop integrals 
and to the various constraints on the deltas. 



15 



7 Results 

The results of the calculations of sections §-^ are presented in figs. |^§ and in tables 
While the gluino sector of the theory is essentially determined by the knowledge of the 
gluino mass (i.e. Mgi), the chargino one needs two more parameters (i.e. M2, /i and tan/?). 
Moreover it is a general feature of the models we are studying the decoupling of the SUSY 
contributions in the limit of high sparticle masses: we expect the biggest SUSY contributions 
to appear for such masses chosen at the lower bound of the experimentally allowed region. 
On the other hand this considerations suggest us to constrain the three parameters of the 
chargino sector by the requirement of the lighter eigenstate not to have a mass lower than the 



experimental bound of about 70 GeV [g^ . The remaining two dimensional space has yet no 
constraint. For these reasons we scan the chargino parameter space by means of scatter plots 
for which M,g = 250 GeV, M^ = 50 GeV, 80 GeV < |/i| < 300 GeV and 2 < tan/3 < 30; for 
every choice of these two parameters, M2 is determined imposing to the lighter eigenstate a 
mass of about 70 GeV. In the plots we sum all contributions coming from different graphs 
proportional to a common mass insertion (the actual values of the coefficients are obtained 
multiplying the points in the plots by the MI). 

In the tables we report the contribution of each diagram and the explicit dependence on 
the mass insertion parameters. We evaluate the coefficients varying Mgq and Mgi between 
250 GeV and 1 TeV. The other parameters in tab. ^ are fixed from the scatter plots in order 
to give the best SUSY contributions to Cg and Cio- 

Thus, with /i ~ -160, Mgi ~ M,, ^ 250 GeV, Mo ~ 50 GeV, M^^ = 90 GeV, tan/3 ~ 2 
one gets 

/ C^'^^Mb) = -I.2i5-,)LL + Om{S^,)LR - 0.51(5|)z.L .4.. 



In order to numerically compare eq. ( ^q) with the respective SM results we note that the 
minimum value of [Cg^{s)] {Mb) is about 4 while C^q'^ = —4.6. Thus one deduces that 
SM expectations for the observables are enhanced when C^^\Mb) is positive. Moreover the 
big value of Ciq^{Mb) implies that the final total coefficient Ciq{Mb) can have a different 
sign with respect to the SM estimate. As a consequence of this, the sign of asymmetries can 
be the opposite of the one calculated in the SM. 

The diagonal contributions to Cg, Cio introduced in sect.^, and computed in the same 
range of the parameters (with Mh+ chosen just above the experimental threshold of about 
60 GeV) are 

|Cg^-^(M^) = -0.35 , . 

\ctr{MB) = -Q.21 ^ ' 

The sign and the value of the coefficient C^ has a great importance. In fact the integral 
of the BR (see eq. (^) is dominated by the \C'!\'^ /s and CyCg term for low values of s. In 
the SM the interference between O7 and Og is destructive and this behavior can be easily 
modifed in the general class of models we are dealing with. 
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In the following, according to the discussion of sect. P, we give the configurations of the 
various Ss for which we find the best enhancements and suppressions of the SM expectations. 



Best enhancement. 










cr'iMB) 




■-^10 


{^2t)LL 


{5^23)lR 


R 


0.41 


1.5 


-8.3 


-0.5 


0.9 


Afb 


0.41 


0.96 


-2.1 


-0.5 


0.15 


Afb 


0.28 


0.96 


-2.1 


-0.5 


0.15 



It is important to note that with such choices the behavior of the asymmetries in the 
low s region of the spectrum is greatly modified: the coefficients of the operators Q-j 
and Qg sum up instead of cancel each other in such a way that the asymmetries are 
never negative. It is also important to stress that the asymmetries get their extremal 
value with a rather small {523)lr- the enhancement given here will survive possible 
future constraints on this insertion. 



Best enhancement with C7 


<0. 








C'r^'iMs) 






('^23'*) LL 


{^^23)lR 


R, 


-0.41 


1.5 


-8.3 


-0.5 


0.9 


Afb, Afb 


-0.28 


0.75 


0.36 


-0.5 


-0.15 


Best depression. 




Cr'iMB) 


riMi 


r'Mi 


{^23) LL 


{5l3)LR 


R 


-0.28 


-1.3 


5.8 


0.5 


-0.6 


Afb-, Afb 


0.28 


-1.5 


8.3 


0.5 


-0.9 



The plots of BR{s), Afb{s) and Afb{s) are drawn in figs.|[||. Here both SM and SUSY 
results are shown. The discontinuity in the Afb plot at s = 0.7 corresponds to the point at 
which we have stopped the corrections 0{l/ml). In fact a model independent description 
of the differential asymmetry in the region 0.7 < s < 0.93 beyond the parton model is still 
lacking. Further the peak wich occur at s = {2mc/mh)'^ ^ 0.3 is due to the perturbative 
remnant of the cc resonance. 

The integrated BRs and asymmetries for the decays B — >■ Xse~^e~ and B — > XgH'^ iJi~ in 
the SM case and in the SUSY one (with the above choices of the parameters) are summarized 
in tab.^. There we computed the total perturbative contributions neglecting the resonances; 
these occur in the intermediate range of the spectrum {J/ip at 3.1 GeV (s = 0.42) and ip' 
at 3.7 GeV (s = 0.59) plus others at higher energies). However it is possible to exclude the 
resonant regions from the experimental analysis by opportune cuts and to correct the effects 
of their tails in the remaining part of the spectrum. 

The results of tab. Q must be compared with the experimental best limit which reads 



BRe^p < 5.8 10 



-5 



(48) 



A comment on the CMSSM (Constrained MSSM) prediction for the observables we have 

g. In this 



computed is now necessary. An analysis on the subject is presented in ref. 
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Observable 


SM 


SUSY 
maximal 


SUSY/ 
SM 


SUSY 
minimal 


SUSY/ 
SM 


SUSY 
(C, < 0) 


SUSY/ 
SM 


BR{e) 


9.6 10"^ 


4.4 10"^ 


4.6 


3.9 10-6 


0.41 


3.9 10-5 


4.0 


Apsie) 


0.23 


0.33 


1.5 


-0.18 


-0.78 


0.31 


1.4 


ApBie) 


0.071 


0.24 


3.3 


-0.19 


-2.7 


0.11 


1.5 


BRifi) 


6.3 10-6 


4.0 10-5 


6.3 


1.610-6 


0.26 


3.4 10-5 


5.4 


^FB(/i) 


0.23 


0.33 


1.5 


-0.18 


-0.78 


0.31 


1.4 


Apsifi) 


0.11 


0.27 


2.5 


-0.27 


-2.4 


0.15 


1.3 



Table 2: Integrated BR, Aps (ind ApB in the SM and in a general SUSY extension of the 
SM for the decays B -^ X^e^e- and B —>■ Xs^^^^. The second and third columns are 



■ieff 



■ieff 



the extremal values we obtain with a positive C/ while the fourth one is the Cj^-' < 
case. The actual numerical inputs for the various coefficients can he found in the text. The 
BR is just the integral of R{s) multiplied by the BR of the semileptonic dominant B decay 
(BR{B -^ X^ep) = 0.105;. 



paper the authors show that the effect of CMSSM on the integrated BRs, considering only 
contributions to Cg and Cio, varies between a depression up to 10% and an enhancement 
of few percents relative to the corresponding SM values. The asymmetries get even smaller 

w) yields in 



corrections. On the other hand a direct computation of Cj^^^^^{M\ 



- 0.59 < C 



MSSM 

7 

MSSM/ 



-0.26<Cf^'""(M, 



(Mw) < +0.49 in the large tan/3 regime, 

w) < —0.20 in the low tan/3 regime. 



(49) 



It is worth noting that comparing the above intervals with the experimentally allowed region 
obtained via RG evolution at the Mw scale of the limits in eqn. (53) (we use only the SM 
contribution to Cg', the inclusion of the MSSM corrections does not change significantly the 
result) 

- 0.39 < CriMw) < -0.099 and 0.66 < CjiMw) < 0.95 (50) 

it is excluded that the CMSSM could drive a positive value for Cj {Mb). For what concerns 
the negative interval of values of Cj {Mb) we see that it can be accommodated both in the 
CMSSM and in our framework. 



Looking at figs ^[|| and table |^ we see that the differences between SM and SUSY pre- 
dictions can be remarkable. Moreover a sufficiently precise measure of BRs, Aps^ and Aps^ 
can either discriminate between the CMSSM and more general SUSY models or give new 
constraints on mass insertions. Both these kind of informations can be very useful for model 
building. 



Diagram 


Msq 


M,i 


Cr 


c. 




250 


250 


~Q.l92{5i^)LL-'i^A{5i^)LR 


-0.5l3iS!^,)LL 




250 


500 


-0.125(4)LL-31.2(52^3)iR 


-0.189{6^,,)ll 


^7 — lins 


500 


500 


-0.0449(<5f3)LL - 15.6(5|)lh 


-0.l2{6i,)LL 




250 


1000 


-0.0344(<52^3)z.L - l^.?>{5i^)LR 


-0M63{6^,,)ll 




500 


1000 


-Qm9l{5t^)LL-l^H5i^)LR 


-0M39{6I,)ll 




1000 


1000 


-0.0105(<5f3)LL - 1.2Q{5t^)LR 


-0.0279iSi)LL 








C,o 


Cg 




250 


250 


-10.2{6',,)USi,)RL 


0.763{6i)LR{Si,)RL 




250 


500 


-17. 3{6i,)LR{Si,)RL 


1.29{6i)LR{Si,)RL 


gZ — 2ins 


500 


500 


-9A9{6',,)LRiSi,)RL 


0.712{6i,)LR{Sl,)RL 




250 


1000 


-17.6{5:j,)U5i3)RL 


1.32(5^,,) USi3)RL 




500 


1000 


-16.1{6l,)LR{Si,)RL 


1.21{6i)LR{Si,)RL 




1000 


1000 


-8.85{S',,) LRiSis)RL 


0M4{6',,)lr{SI,)rl 



Table 3: Contributions to the coefficients Cj, Cg and Cio from diagrams involving gluino 
loops. Mgi and Mgg both vary between 250 GeV and 1000 GeV. Exchanging L with R in the 
mass insertions we get the contributions of gluino diagrams to C^, Cg and C[q. For further 
explanations see the caption in tab. H. 
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Diagram 


Msq 


Cj 


c. 




250 


0-MS23)lL 


1-4(<52"3)ll 


WW-f - lins 


500 


0.12{6^^)ll 


0.76((52"3)ll 




1000 


0.033(52"3)ll 


0.32((5«3)ll 




250 


-2.1((52"3)LL-0.25((52"3)LiJ 


-^■n{5^3)LR 


HW-f - lins 


500 


-l.m,)LL-0.27i6^,)LR 


-0.87((5«3)lk 




1000 


-0A5{61^,)ll - 0.27{S^,)lr 


-Qm{5i,)LR 






Cio 


c. 




250 


l-m3)LR{S^3)RL 


-Qm2{5^,,)LR{5l,)nL 


WWZ - 2ins 


500 


mS^3)LR{5i;s)RL 


-0.12(5«3)LH(53"3)^i 




1000 


2.m,)LR{S^,)nL 


-Q.U{5^,)lr{51,)rl 




250 


-8.4((5«3)lh 


0.56((5«3)Lfi 


WHZ - lins 


500 


-IHS^3)lR 


0.74(<52"3)i^ 




1000 


-IHS23)lR 


0.84(52"3)lr 




250 


-0.91(5«3)ll 


0.06((52"3)ll 


WWZ - lins 


500 


-0.47((52«3)ll 


0.031((52"3)ll 




1000 


-0.19{61^,)ll 


0.013(52"3)ll 




250 


2.7(52«3)ll 


-2.7(<52«3)ll 


boxW^ — lins 


500 


1-3(52"3)ll 


-1-3(52"3)ll 




1000 


0.55{6^,)ll 


-0.55(52"3)ll 




250 


-0.97i6^^)LR 


0.97(52"3)LiJ 


hoxHW - lins 


500 


-i-m3)LR 


1-W23)lR 




1000 


-i-m3)LR 


l-2(<52"3)Li? 



Table 4: Contributions to the coefficients C7, Cg and Cio from diagrams involving chargino 
loops. We assume /i = -160 GeV, M2 = 50 GeV, tan/5 = 2, m^ = 50 GeV, M^^ = 90 
GeV while Mgq varies between 250 GeV and 1000 GeV. In the ffist column we indicate the 
number of mass insertions present in each squark line, which charginos are present at the 
vertexes and the kind of graph computed (7-penguin, Z-penguin or box diagram). 



20 




300 




300 



Figure 3: (^23)^^ (above) and [623)11 (below) contributions to Cg coming from chargino 
diagrams as functions of fi (expressed in GeV). Msq is fixed to 250 GeV while tan/3 varies 
between 2 and 30. 
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300 




300 



Figure 4: {52^lr (above) and (^23)^^ (below) contributions to Cio coming from chargino 
diagrams as functions of fi (expressed in GeV). Msq is fixed to 250 GeV while tan/3 varies 
between 2 and 30. 
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Figure 5: Differential branching ratio for the decay B — >■ X^^^^ . The solid line corresponds 
to the SM expectation; the dashed and dotted-dashed lines correspond respectively to the 
SUSY best enhancement (Cj^^ = 0A1,C^" = 1.5, CfJ^ = -8.3; and depression (Cf^ = 

5.8); the dotted line is the maximum enhancement obtained 



-0.28,C^ 



MI 



-1.3, C 



MI 
10 



Without changing the sign of Cj (C^^^ = -0.41,Cg^^^ = 1.5, C*^^ = -8.3;. 
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Figure 6: Forward-Backward asymmetry (Aps) for the decay B -^ Xst^l . The solid line 
corresponds to the SM expectation; the dashed and dotted-dashed line corresponds to the 
SUSY best enhancement (Cj^^ = OAl,Ci" = 0.96, C^^^^ = -2.i; and depression (C^^^ = 
.28,Cq^^ = —1.5,CiQ^ = 8.3j; the dotted line is the maximum enhancement obtained without 



changing the sign of Cj (C^ 



0.75, C,^^ 



0.36;. 
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Figure 7: Forward-Backward asymmetry (Aps) for the decay B -^ Xsfi~^fi . The solid line 
corresponds to the SM expectation; the dashed and dotted-dashed lines correspond respec- 

= 0.28,C^^^ = 0.96, Ci^^ = -2.i; and depres- 
the dotted line is the maximum enhancement 



tively to the SUSY best enhancement (Cj 



sion (Cj 



0.2S,C, 



MI 
9 



-1.5, C^i^' 



U); 



obtained without changing the sign of Cy (C^^^ = -0.28,0^" = 0.75, Cf^^ = 0.36;. 
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Figure 8: Forward-Backward asymmetry (Aps) for the decay B^ -^ X^e'^e (See the caption 
in fig. 0j. 
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8 Conclusions 

In this paper an extensive discussion about SUSY contributions to semileptonic decays B -^ 
Xse^e~, B -^ XsyU+yU^ is provided. We see that the interplay between b -^ s'-f and B -^ 
Xsi^i" is fundamental in order to give an estimate of the SUSY relevance in these decays. 
The two kinds of decays are in fact strongly correlated. 

Given the constraints coming from the recent measure of 6 — > 57 and estimating all 
possible SUSY effects in the MIA framework we see that SUSY has a chance to strongly 
enhance or depress semileptonic charmless B-decays. The expected direct measure will give 
very interesting informations about the SM and its possible extensions. 
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